AD-A260  470 


MENTATION  PAGE 


form  Approved 
0MB  No  0704-0188 


.  to  '  ♦'our  otr  <nc'ua*nq  ttme  io^  fev*ev%inq  inMfuaic''%  vearcninq  aata  soufc«, 

q  and  fewte'^»r»q  the  co^eaion  of  «nfof*wation  Sena  commenii  reqa^ding  tni\  Durden  e^tlmate  0^  an»  other  aspect  O*  tn-s 
iq  this  ouraen  to  wVashrn<jton  »^eadduaners  Serv/ces.  Directorate  w'  )nfo''mat>on  Operations  and  >‘*‘tx5rts.  12  tS  ;eMerson 
to  the  O^'ce  of  Management  and  Budget  Paperwork  Heducticn  Protect  (0  ?04-0 188)  Aashington  DC  20S03 


1.  AGENCY  USE  ONLY  (Leave  blank)  I  2.  REPORT  DATE 


4.  TITLE  AND  SUBTITLE 

'^'2.  C^'»  V  C  V'  tvW  I  A  IX-oi 
LiaDl^ 


3.  REPORT  TYPE  AND  OATES  COVERED 


S.  FUNDING  NUMBERS 


6.  AUTHOR{S) 


Ho  tjL(AA  K  .  i-<Ki. 


7.  PERFORMING  ORGANIZATION  NAM£(S)  AND  AOORESS(ES) 

~vC^  e.  C-\  Clc  ( ( (2. «  A  ^ — 1 

'  U  /o.U  ^  \ 

H  V  cD_  \Cc5v 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


00(>i3 

OTIC 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  In  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
osltlon,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 


Abstract 


93-03406 


Inversion  of  light  scattering  by  a  .single  spherical  scattercr  with  or  without  a  layered 
structure  is  studied.  By  numerical  simulation  of  the  statistical  properties  of  experiments, 
we  can  select  those  angles  that  will  maximize  the  ability  of  the  experiment  to  resolve  a 
sphere  from  a  layered  scattcrer.  Applying  statistical  decision  theory,  we  obtain  a  criterion 
for  decision  making  that  minimizes  the  probability  of  incorrect  guesses. 


14.  SUBJECT  TERMS 

VMSvse. 


oiA 

\y\  Jiv 


IS.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  118.  SECURITY  CLASSIFICATION  119.  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 

OF  REPORT  I  OF  THIS  PAGE  I  OF  ABSTRACT  I 


unclassified 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298  (Rev  2-89) 

PrescfibW  by  ANSI  Sid  ^39-18 
298  102 


nr/. 


Acce^io  1  ►O'  _ ^ 

MIS  CRA&l  iS 

DliC  tab  □ 

-  '.(.-d  U 


QUASI-BINARY  DECISION  MAKING  USING  LIGHT  SCATTERING 


Po  Hu  and  M.  Lax 

Physics  Department,  City  College  of  New  York 
New  York,  New  York'^lOOSl 


^  jii  iLj’nTy  '  1  dts 

t  A,/  111  .1  ‘d  or 
"^t  SufC'-jl 


RECENT  PUBLICATIONS,  SUBMITTAL  FOR  PUBLICATION  AND 
PRESENTATIONS: 


’P.  Hu,  M.  Lax,  “Single  Scattering  Inversion  Problem’’,  in  Proceedings  of  the  1990 
CRDEC  Scientific  Conference  on  Obscuration  and  Aerosol  Research. 


^B.  Yudanin,  M,  Lax,  “POST  Adaptation  for  a  Numerical  Solution  of  the  Spherically- 
Symmetric  Riemann  Problem,’’  J.  Mod.  Phys.  285-298  (1990). 

^B.  Yudanin,  P.  Hu,  M.  Lax,  “Numerical  Solution  of  the  Riemann  Problem  in  the 
Presence  of  an  External  Energy  Source”,  in  Proceedings  of  the  1988  CRDEC  Scientific 
Conference  on  Obscuration  and  Aerosol  Research. 


'^M.  Lax,  B.  Yudanin,  “Early-Time  Hydrodynamic  Response  to  Pulsed  Laser  Radiation,” 
Digest  -  Inti.  Conf.  on  Optical  and  Millimeter  Wave  Propagation  and  Scattering  in  the 
Atmosphere,  Florence,  Italy,  1986. 

^B.  Yudanin,  M.  Lax,  “Hydrodynamical  Re,sponse  to  Uniform  Laser  Absorption  in  a 
Droplet,”  in  Proceedings  of  the  CRDC  1985  Scientific  Conference  on  Obscuration  and 
Aerosol  Research,  edited  by  R.  H.  Kohl  and  D.  Stroud  (1985). 


Abstract 

Inversion  of  light  scattering  by  a  single  spherical  scatterer  with  or  without  a  layered 
structure  is  studied.  By  numerical  simulation  of  the  statistical  properties  of  experiments, 
we  can  select  those  angles  that  will  maximize  the  ability  of  the  experiment  to  resolve  a 
sphere  from  a  layered  scatterer.  Applying  statistical  decision  theory,  we  obtain  a  criterion 
for  decision  making  that  minimizes  the  probability  of  incorrect  guesses. 
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Inlrodiictioii 


Inversion  prohlenis  have  existed  in  various  branches  of  engineering  and  physics  for 
a  long  time,  but  in  the  past  twenty  years  they  have  receiverl  far  more  attention  than  ever 
before  because  of  tlie  available  of  high  speed  computers.  !n  the  present  paper,  we  are 
going  to  invert  light  scattering  information  to  recognize  an  inner  structure  of  a  spherical 
object.  Because  intensities  of  scattering  light  are  highly  nonlinear  functions  of  the  size  and 
index  of  the  scatterer,  and  noise  is  present,  arbitrary  pattern  recognization  is  difficult.  We 
will  restrict  ourselves  to  distinguishing  two  kinds  of  refraction  index  patterns  (uniform  or  a 
layered  sjdiere)  with  a  range  of  parameters.  We  assume  that  the  scatterer  may  be  one  of 
two  kinds  of  objects;  a  sphere  with  a  uniform  refractive  index  or  a  layered  sphere  with 
different  refraction  indices  for  the  core  and  shell,  re.spectively.  Our  problems  arc:  (1)  For 
given  experimental  conditions,  is  there  enough  information  to  make  a  decision?  What  is 
the  best  choice  of  angles  to  yield  the  most  significant  statistics.  (2)  For  given  experimental 
data,  how  should  one  make  a  decision  that  minimizes  wrong  guesses. 

We  assume  that  we  already  know  the  following  facts:  The  wave  length  X  of  the 
light  in  vacuum  is  .4416  microns.  For  the  u=  uniform  sphere  (hypothe.ses  1),  the  varying 
parameters  are  the  radius  R“  and  refraction  index  in  the  ranges  4</?“/X<8  and 
1.33</t“  <  1.8.  For  the  /=  layered  .sphere  (hypothe.ses  2),  the  varying  parameters  are  the 
inner  radius  /?(„,  core  and  shell  refraction  indices  /t|n  t'titl  They  are  in  the  ranges 

4<r{„/X<1,  1.33<n|n^l.5  and  1..55 <«m,i  <  1.8.  The  outer  radius  /?oiit  of  the  layered 
sphere  is  fi.xed  at  Moreover,  we  assume  there  is  a  Gaussian  noise  added  to  the 

scattered  intensity  with  a//=0. 1,  where  o  is  the  width  of  the  intensity  distribution  of  the 
Gaussian  noi.se  and  /  is  the  mean  intensity  of  the  light  scattered  at  a  given  angle. 

The  method  of  examining  general  experimental  data  to  decide  between  two 
hypotheses  is  a  classic  problem  in  decision  theory.  A  procedure  in  the  absence  of  a  priori 
information  was  propo.sed  by  Neyman  and  Pearson’  in  1933.  Their  results  are  expres.sed 
in  terms  of  “maximum  likelihood  ratios”.  An  excellent  overview  is  given  by  Kendall  and 
Stuart^.  A  readable  description  of  the  Bayes  theorem  approach  to  the  same  problem  is 
given  by  van  Trees^.  The  close  connection  between  these  two  approaches  is  touched  on 
by  Middleton"*  in  his  section  on  binary  detection  systems. 

Our  problem  is  more  general,  in  that  we  must  estimate  some  continuous  parameters 
first,  in  order  to  make  the  best  binary  decision.  In  the  present  application,  these  parameters 
are  radii  and  indices  of  refraction.  Thus  our  problem  mixes  continuous  parameter 
estimation  with  discrete  parameter  detection. 

Distinguishability 

We  propo.se  to  study  the  statistical  properties  so  that  we  may  obtain  a  crit  -rion  for 

no 


measuring  tlio  distingiiishabiliiy. 


Suppose  the  experimental  light  scattering  data  have  been  taken  at  (spherical) 
scattering  angles  0,  and  tj)=();  the  intensity  of  the  scattering  light  per  unit  solid  angle  is 
denoted  as  /(0,).  We  want  to  relate  the  observation  infomiation  to  the  structure  of  the 
scatterer.  Because  noise  is  always  involved  and  can  not  be  separated  from  the 
experimental  data,  it  will  mix  the  two  hypotheses  such  that  they  can  not  be  distinguished 
when  the  noise  level  is  relatively  high.  The  noise  level  depends  on  the  experimental 
apparatus,  the  environment  and  the  kind  of  data  taken.  For  a  given  apparatus  and 
environment,  we  should  measure  such  data  to  minimize  the  relative  noi.se  level. 

For  a  single  scattering  process,  the  data  measured  at  different  angles  are  not 
independent;  they  are  correlated  through  complex  formulas  (Mie  scattering  for  a  sphere  and 
shell).  Because  of  the  complex  relationships  of  the  scattering  formulas,  we  can  not 
analytically  solve  and  will  use  numerical  simulations  to  our  problem. 


Let’s  say  the  object  is  a  uniform  sphere  with  the  parameter  element  r\^^^e{R‘‘ } 
in  the  allowed  range  and  the  experimental  obserx'ation  is  taken  at  M  angles  Bj,  .  .  .  ,0\/. 
Here  A/  is  a  moderate  number  of  order  10.  We  can  regard  the  scattering  as  a  mapping 
from  the  parameter  space  to  an  M  dimensional  observation  space 

L'''^={/(0i),....,/(0.v/)}.  Because  noise  is  present,  one  point  in  rj  space  maps  to  an  M 
dimensional  "box"  in  space.  A  similar  mapping  also  applies  for  a  layered  object.  If,  at 
each  angle  0,-,  separate  measurements  with  P  different  polarization  intensities  are  made, 
corresponding  to  different  polarization  of  the  incident  and  scattered  beam,  then  PM  will  be 
the  dimension  for  the  observation  space.  For  simplicity  of  notation,  however,  we  do  not 
always  introduce  an  explicit  polarization  index. 

In  the  following,  we  will  define  a  resolution  criterion.  First,  we  generate  a  set 
{ri^^^}  which  includes  N  elements  of  random  parameters  in  the  allowed  ranges  for 
hypothesis  j  (J=\  for  uniform  sphere  and  j=2  for  layered  object),  where  m=l,...A/.  Here  N 
is  a  large  number  of  order  1000.  The  m\h  element  rim^  of  uniform  .sphere  set 
has  two  parameters  and  and  the  mth  element  nlp  of  the  layered  object  set 
contains  three  parameters  «L,in  and  With  the  added  noise,  one  image  for 

each  of  the  two  sets  and  are  obtained  in  the  observation  space  Within 

the  overlap  region  of  the  two  images  in  the  observation  space  the  two  kinds  of  objects 
are  indistinguishable  from  the  given  M  intensity  measurements.  We  will  give  the  definition 
of  overlap  later.  An  event  producing  an  image  in  the  overlap  region  of  the  observ’ation 
space  is  regarded  as  an  indistinguishable  event  for  tlie  two  hypotheses.  Counting  the 
number  of  the  events  which  overlap  in  the  observation  space,  we  may  get  a  measure  of 
indistinguishability  between  the  two  cases.  The  ratio  of  the  number  of  overlapped  events  to 
the  total  number  events  repre.sents  a  measure  of  the  indistinguishability. 

We  study  the  probability  properties  of  the  two  hypotheses.  In  the  observation  space, 

m 


an  imago  point  from  one  hypotlicsis  can  always  have  a  probability  of  overlap  with  the  an 
image  of  the  other  kind  of  hypothesis.  Therefore  we  need  to  precisely  deline  overlap  of 
images  in  the  observation  space.  We  use  N(l{Q))  for  the  number  of  events  such  the 
intensity  is  in  the  range  of  /(0)  to  /(OHJ/fO)  and  N  for  the  total  number  of  events.  We 
then  define  an  overlap  function: 


I 


N 


N 


(1) 


where  the  subscript  p  denotes  the  polarization  and  the  super.script  (1)  and  (2)  are  for 
different  hypotheses.  The  set  of /V  (-lOOO)  points  in  po.s.sesses  a  subset 

which  overlap  in  the  sense  that  it  may  no  longer  be  possible  to  distinguish  whether  case 
j=\  or  7=2  is  the  correct  with  respect  to  a  single  measurement.  We  can  choose  a 
threshold  a,  ^  of  indistinguishability  for  0,  and  measured  intensity  Ip  by  using  the  criteria 


tor  all 

to  determine  the  elements  of  subset  [ri*'t(Q^)|  por  obserc'ations  at  many  angles,  the  joint 
set  of  subsets  for  all  measured  ancles 


,(-)  I,' 


(2) 


[Ti0)(0,,...,0„)l=nlri'‘\0,)l 


(1), 


(3) 


measures  the  overall  fuzziness  of  the  experiment.  As  an  application,  we  assume  the 
measurement  is  taken  at  given  angles  0=25,40.90,  105,  125,140  for  two  polarizations  (12 
measurements).  For  simplicity  we  set  all  a,-,^  to  be  the  same  a,  p=a.  In  Fig.  1  we  show 
the  original  set  (dots)  of  parameters  for  the  unifomi  sphere  chosen  by  Monte  Carlo 
techniques,  and  its  joint  subset  for  12  measurements  (triangles).  The  ratio  of  the  numbers 
of  the  elements  of  the  joint  sub.set  and  the  original  set  is  about  10%.  Therefore,  about 
90%  of  uniform  sphere  events  (the  dots  not  covered  by  the  triangle  in  Fig.  1)  are 
distinguishable  from  the  layered  events  in  the  space.  The  superscript  12  is  the  number 
of  dimensions  of  the  observation  space,  that  is  the  number  of  measurements.  Tlie 
remaining  10%  of  the  events  for  the  uniform  sphere  (the  triangles  in  Fig.  1)  are 
indistinguishable  from  the  layered  scatterer  for  the  given  set  of  12  measurements. 

To  produces  a  better  resolution  between  these  two  hypotheses,  a  trivial  approach  is 
to  increa.se  the  number  of  the  detectors  so  that  the  events  in  the  new  space  (A/>12) 
will  not  overlap  as  much  as  in  the  original  space  In  most  experiments, 

distingui-shability  is  limited  by  the  number  of  detectors.  To  make  the  most  of  the 
equipment  for  a  better  resolution,  one  can  rearrange  the  detectors  for  some  optimized 
angles  so  that  the  number  of  the  elements  in  the  joint  subset  is  minimized.  This  can  be 
done  by  the  following  iteration  procedure:  First  we  randomly  generate  two  parameter  sets 
{T|0^}  and  as  initial  sampling  sets.  Use  Eq.  (4)  below  to  find  the  best  angle  0i  by 

taking  a  minimum  of  the  overlap  function  F  for  two  polarizations  over  Oj .  Substitute  ©i 
and  the  original  sampling  sets  into  Eq.  (5),  below,  to  obtain  the  subsets  [n^'^(8i)]  and 
[Ti^^^(0i)]  by  Eq.  (5).  Iterating  this  procedure  by  u.sing  the  new  sets  [ti^0(0i)]  in  Eqs.  (4) 
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and  (5').  we  will  obtain  the  second  best  angle  0i  and  the  subset  ,02)1.  Repeat  the 

iteration  procedure  (4)  and  (5)  for  angles  0|,02,0.w  •  •  •  until  the  number  of  elements  of 

the  joint  subset  |ti^‘'(0i . ,0,i)|  is  less  than  a  desired  value.  With  the  optimization 

completed  for  i-1  angles,  we  can  optimaize  over  0,  using: 

->  oa 

min  n  . ,0,_, )|,[n‘“’(e, . 0,_i)l)  (4) 

p  =  \  0 

. ,e,-t.e,)i  if  F^(0„/^;(Ti‘'\0i,....,0,-i)Mn‘‘'(0i . ,0,-1 

for  all  (5) 

Here  we  assumed  that  intensities  of  two  polarizations  (the  subscript  p)  for  each  angle  have 
been  measured. 

Decision 

In  this  section,  we  apply  statistical  decision  theory  to  the  inversion  problem.  Let  // 1 
and  H  2  denote  the  two  hypotheses;  1  for  the  unifonn  sphere  and  2  for  the  layered  sphere. 
Suppo.se  we  have  obtained  a  set  of  experimental  data  for  the  intensities  at  several  angles 
fcxp(0<)-  We  want  to  decide  which  class  (uniform  sphere  or  layered  object)  the  scatterer 
belongs  to.  We  can  use  the  least  square  fit  to  find  the  best  fit  for  the  hypothesis  h: 

^/:  =  'jjnl(/cxp(0,)-/(eoTl^'’^)]^  (6) 

within  the  pemiitted  parameter  space.  A  simple  statistic  to  decide  between  these  cases  can 
be  chosen  as  V1-V2.  Suppose  the  experimental  data  are  from  a  uniform  sphere,  /:=!. 
Then  vj  is  dominated  by  the  experimental  noise,  usually  a  small  value,  while  V2  is  a  large 
value  because  it  is  not  dominated  by  noise  but  by  the  shift  because  of  an  incorrectly 
chosen  hypothesis.  For  some  experimental  data  the  fitting  with  wrong  hypothesis  may  be 
small  if  the  number  of  the  detectors  is  not  large  enough.  The  wrong  fitting  value  of  v 

covers  a  large  range  of  a  uniform  scale  space.  Therefore  we  use  a  log  scale  for  the  statistic 

/?=logio(Vi/V2) 

to  decide  between  the  two  hypotheses.  The  decision  rule  can  be  obtained  as  follows. 
Let’s  assume  we  know  the  the  conditional  probability  P{R  1///,)  of  getting  R  under 
hypothesis  h.  According  to  Neyman-Pearson’"'*,  when  a  priori  probability  and  the  cost  of 
the  decision  are  unknown,  we  may  use  a  constraint  condition  on  the  probability  Pp,  of  a 
"false  alarm"  is: 

Pp=a^j°°PiR  iH^dR  (7) 

to  find  the  threshold  X,  where  a  is  the  value  permitted  for  a  false  alarm  (we  say  H2  while 
H\  is  true).  After  finding  the  threshold,  we  shall  make  decision  by  the  criterion: 
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if  \{R)>X 

choose  Ih 

(<S) 

if  A(R)<X 

choose  // 1 

(9) 

where 

A{R)=P{R  1 

\Hy)/P(R  I//,) 

(10) 

is  the  maximum 

likelihood  ratio. 

The  conditional  probabilities,  P{R  !//>,),  can  be  obtained  (ahead  of  time)  by 
numerical  simulation.  As  an  application,  we  generate  a  set  consisting  of  lOOO  elements  of 
random  parameters  for  each  hypothesis.  We  then  calculate  the  intensities  of  the  light 
scattering  for  these  parameters  at  the  angles  0=90,  105,  ,120,  135,  150,  165  for  both  parallel 
and  perpendicular  polarizations.  Finally  we  add  10%  noise  to  the  calculated  intensities  and 
regard  the  re,sult  as  pseudo  experimental  data.  To  get  the  distribution  protiles  for  both 
cases,  we  also  use  the  least  square  method 


to  fit  the  pseudo  experimental  data  with  the  best  parameter  Here  /’c^p  was  computed 
for  source  i'  (pseudo  experimental  data),  is  the  best  least  square  fit  for  a  set  of 
experimental  data  of  source  .v  by  the  hypothesis  It  within  the  permitted  paramete. ..  Becau.se 
we  used  wide  ranges  for  the  parameters,  the  intensities  have  hundreds  of  oscillations  over 
the  varying  parameters.  It  is  difficult  to  locate  a  global  minimum  for  v^,  because  it  has 
hundreds  of  oscillations  over  the  varying  parameters.  To  make  the  programs  more 
efficient,  we  made  lookup  tables  for  the  Bessel  and  Legendre  functions.  Defining 
y?^=logio(Vi  /vj)  and  counting  the  number  of  the  events  in  which  /?  '  falls  into  the  interval 
(R,R+dR),  denoted  as  /V(/?'),  we  obtain  the  distribution  profile  for  the  source  s.  In  Fig.  2, 
N(/?^)  vs  R^  is  plotted  for  both  sources.  The  left  profile  of  this  figure  is  for  the  source 
5=1  the  unifonn  sphere,  while  the  right  one  is  for  the  source  .v=2  the  layered  object.  The 
small  overlap  between  the  two  curves  in  Fig.  2.  shows  that  the  resolution  of  these  two 
cases  is  quite  good. 
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Fig.  1.  Joint  sub.set  and  original  parameter  set  for  the  uniform  sphere.  The  vertical  scale 
is  for  R‘*/X.  The  1000  dots  (and  triangles)  are  randomly  .selected  from  the  available 
parameter  set  for  the  uniform  sphere.  The  triangles  i\rc  for  the  elements  of  the 

joint  subset  of  uniform  sphere  for  the  6  given  angles  and  two  polarizations;  they  are  the 
points  indistinguishable  from  the  layered  objects.  The  ratio  of  the  numbers  of  the 
triangles  and  dots  is  about  10%,  or  about  90%  of  the  uniform  sphere  events  are 
distinguishable  from  layered  objects.  The  original  parameter  set  consists  1000  elements, 
each  one  has  two  components  of  random  numbers  for  the  parameters  of  the  uniform 
sphere.  The  two  parameters  are  the  radius  and  refraction  index  in  the  ranges  of 
4^“/X<=8  and  1.33<«“<1.8.  The  measurements  are  taken  for  the  intensities  of  two 
polarizations  (parallel  and  perpendicular)  at  6  angles  0=25,  40,  90,  105,  125,  140. 


115 


100 


N(R') 

50 


0 

-5  0  5 

/?'=log,o(vl /V2) 

Fig.  2.  Distribution  profiles.  The  vertical  axis  is  /V(/?  ')  the  number  of  events  per  unit  R, 
dR=0A  is  used.  The  left  profile  is  the  distribution  profile  for  the  appropriate  for  uniform 
spherical  source,  while  the  right  one  is  for  the  layered  scatterer.  The  measurements  are 
taken  for  the  intensities  of  two  polarizations  (parallel  and  perpendicular)  at  6  angles 
0=90,  105,  120,  135,  150,  165  with  10%  noise.  The  uncertainty  of  the  parameter  ranges 
are  follows:  For  the  sphere,  the  two  parameters  arc  the  radius  and  refraction  index  in  the 
ranges  of  4^“/A.<8  and  1.33<«“<1.8.  For  the  layered  object,  the  three  parameters  are 
the  inner  radius  and  refraction  indices  in  the  ranges  of  1.33</i'n<1.5  and 

1.55<n-n^l.8. 
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